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Résumé: Cette note introduit une classe de matrices dont les déterminants 
sont faciles à calculer. L 'exemple le plus frappant est obtenu en considérant 
la matrice entière symétrique M(n) avec coefficients Mij E {0, 1}, < i, j < 
n définis en considérant les coefficients hinomiaux modulo 2 et en posant 
Mij = (mod 2). Le déterminant det(M(n)) est alors étroitement 

relié à la suite de Thue-Morse comptant les coefficients non-nuls d'entiers 
binaires. 

1 Définitions et résultats 

Dans ce papier, la lettre b dénotera toujours un entier > 2 et la notation 
n = ^fj6*, < < 6 — 1 désignera un entier non- négatif n écrit en base b. 
Définition 1.1. Une matrice (finie ou infinie) 

M{n) = (M M ) Q < M . <n , n G N U {oo} 

(à coefficients Mjj dans le corps ou anneau unifère commutatif favori du 
lecteur) est b— autosimilaire si on a M^o = 1 et 



pour tous < s = J2 o~ib l , t = J2 T ib l < n. 

Une matrice M(n) qui est b— autosimilaire pour b < n est entièrement 
déterminée par la sous-matrice M = M(b) de coefficients Mij = Mij, < 
i,j < b. On appelera M la matrice de définition de M(n). 

Exemple 1.2. Pour b = p un nombre premier, la matrice symétrique S 
et la matrice triangulaire inférieure T 

S=(s i , j =^ + . J ^J (modp)j et T=(t m =Q (modp)j 

à coefficients dans le corps fini Z/pZ sont p— autosimilaires (voir section 2). 

Remarque 1.3. Une autre raison de s'intéresser aux matrices b— autosimilaires 
provient du produit tensoriel: Soit (-Bij)o<ij<b une matrice à coefficients 



i 



dans un corps K définissant un endomorphisme de V = K b . La matrice 
B(b d ) (qui n'est b— autosimilaire que pour .60,0 = 1) définie par 

d— 1 d-1 d-l 

B a>t = [] B ai , Ti , 0<s = J2 °ib\ t = nb l < b d 

i=0 i=0 i=0 

définit alors un endomorphisme diagonal 

_g<g> d . y® d > y® d 

vq ® ■ ■ ■ <g) Vi <8> • • • <8> Vd-i i — > Bvq <8> ■ ■ ■ <8> Bvi <8> • • • (g) Bvd-i 

où uo, . . . , «d-l € V = K b . La normalisation i?o,o = 1 permet alors de 
s'affranchir des puissances entières b d de b et définit un analogue d'une 
telle matrice pour toute dimension entière ainsi que pour une dimension 
oo dénombrable. 

Une matrice b— autosimilaire M(n), n G {b, b + 1, . . .} U {oo}, est non- 
dégénérée si les b matrices M(k) = {Mi^Q^j^, k = 1, . . . , b sont toutes 
inversibles. On pose alors d(0) = 1 et 

d(k) = det(M (fc + l))/det(Af (fe)), k = 1 . . . , b - 1 . 

En particulier, une matrice b— autosimilaire triangulaire T(oo) est non- 
dégénérée si et seulement si sa matrice de définition T est inversible. Les 
nombres d(0), . . . , d(b — 1) introduits ci-dessus ne sont alors rien d'autre que 
les coefficients diagonaux de T. 

Théorème 1.4. Soit b > 2 un entier et M = M{n) une matrice 
b— autosimilaire non-dégénérée. 

(i) On a une factorisation unique 

M = LDU 

où L est une matrice b— autosimilaire, unipotente (L^i = 1, < i < n) trian- 
gulaire inférieure, D est b— autosimilaire diagonale et U est b— autosimilaire, 
unipotente triangulaire supérieure. 

(ii) On a À,i = d(i), < i < b et 

n-l 

det(M(n)) = det{D(n)) = {{ ïl d ( K i) + 

fc=0, k=J2^ib i i 

pour tout n G N. 

(iii) L'ensemble des matrices (infinies) triangulaires inférieures (ou 
supérieures) b— autosimilaires et non-dégénérées est un groupe: Si R et T 
sont deux telles matrices définies par R et T, alors RT est définie par RT. 

Remarques 1.5. (i) Une matrice M(n), n G N qui est b— autosimilaire 
et non-dégénérée s'inverse donc relativement facilement: 

(M(n))- 1 = (U(n)T l (D(n)y l (L(n)y l 
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avec U, D, L comme dans l'assertion (i) trois matrices inversibles triangu- 
laires. Leurs inverses sont donc b— autosimilaires et peuvent se calculer à 
partir des inverses U" 1 , D^ 1 , L -1 de leurs matrices de définition. 

(ii) Le théorème 1.4 permet très souvent de calculer det(M(n)) même si 
la matrice b— autosimilaire M = M{n) (avec n € N) n'est pas non-dégénérée: 
Il suffit de considérer la matrice 6— autosimilaire associée à une perturbation 
générique D+tA (avec ^4o,o = 0)) et d'évaluer le résultat en t = 0. L'exemple 
2.4 de la section suivante illustrera ce procédé. 

(iii) Si la matrice B considérée dans la remarque 1.3 est une matrice 
inversible, alors la preuve du théorème montre que la matrice B(b d ) est 
inversible et que son inverse est C{b d ) où BC = id. 

Preuve du théorème 1.4. Soit M la matrice de définition de la matrice 
b— autosimilaire M(n). Comme M[n) est non-dégénérée, on a une factori- 
sation unique (donnée par le procédé d'orthogonalisation de Gram-Schmitt) 

M = LDÛ 

avec L unipotent triangulaire inférieure, D diagonale et U unipotent tri- 
angulaire supérieure. Comme Lo,o = A),o = ^0,0 = 1> les matrices L,D 
et U définissent des matrices b— autosimilaires L(l),D(l) et U(l) pour tout 
l € NU {oo} et un petit calcul montre que ces matrices sont respective- 
ment unipotente triangulaire inférieure, diagonale et unipotente triangulaire 
supérieure. 

On a alors (avec L = L(oo),D = D(oc), U = (oo) et M = M(oo)) pour 

o<s = J2aib\ t = J2nb i 

(LDU) a>t =Y^ L s,kDk,kU k ,t 

k 

k=Y J K i b i * 
6-1 

i Ki=0 

= IiM ai , Ti =M Sit 

i 

ce qui montre l'assertion (i) (modulo l'unicité, laissée au lecteur). 
Les égalités 

det(M (i)) = det(D(0) pour l = l,...,b 

impliquent par récurrence sur < i < b qu'on a l\i = D^i = d(i). On a 
donc 

n— 1 

det(M(n)) = det(D(n) = JJ D Kk 

k=0 
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n— 1 

fe=0, k=Y J K i b i i 

ce qui prouve l'assertion (ii). 

L'assertion (iii) résulte d'un calcul similaire à celui utilisé dans la preuve 
de l'assertion (i), à savoir: 

( RT )s,t = ^2 R s,kTk,t 
k 

k=J2 i 
i Ki=0 



pour R et T deux matrices triangulaires inférieures b— autosimilaires. Il 
faut également vérifier la condition de normalisation (i?T)o,o = 1 qui résulte 
évidemment du fait que R et T sont tous les deux triangulaires inférieures. 
QED 



2 Exemples liés au triangle de Pascal 

Fixons un nombre premier p et considérons la réduction (mod p) du tri- 
angle de Pascal formé des coefficients binomiaux Q) , n, i G N définis par 

(l + *) B = Ç^**€Z[s] • 

L'existence de Pautomorphisme de Frobenius (p : a i — ► aP sur un corps 
de caractéristique p montre qu'on a pour n = J2 ViP % 

(1 + x) n = H{l + x pi ) Vi (modp) 

i 

ce qui implique pour k = J2 K iV l 1 & factorisation 

( E q. i) (ï) =n (*) ( mod ^ 

et prouve que les matrices de l'exemple 1.2 sont bien p— autosimilaires. 
L'identité triviale (1 + x) s+t = (1 + x) s (1 + x) 1 implique l'égalité 
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et fournit la factorisation S = S(n) 



entière S avec coefficients S, 



CV) 



L(n)L(n)* de la matrice symétrique 
j < i,j < n comme produit de la 



matrice unipotente triangulaire inférieure L avec coefficients Lj j = Q) , < 
i,j<n et de sa transposé. La réduction (mod p) de cette factorisation 
illustre l'assertion (i) du Théorème 1.4 (avec D = id et U = LÏ) pour les 
matrices de l'exemple 1.2. 
Considérons la matrice 



A = A(n) 



( 
1 


V : 







définie par 



Aa 



i 




si i = j + 1 , 
sinon 



pour < i,j < n. Un calcul facile montre qu'on a 

L{n) = exp A(n) = JT 



k=0 



(où L(n), donné par Ljj = Q, < i, j < n, est la matrice unipotente 
triangulaire inférieure introduite ci-dessus) ce qui implique que la matrice 
R = L^ 1 = exp (—A) est donné par 



R 



1,3 



(ceci se démontre évidemment aussi aisément en calculant LR). Le calcul 
de la matrice S{n)~ l = i?(n)' R(n),n £ N (dans l'anneau Z des entiers ou 
dans le corps fini Z/pZ) est donc tout à fait explicite. 

Soit p un premier que nous fixons et soit a G Z/pZ un élément du corps 
fini à p éléments. Rappelons que le symbole de Legendre (j^J est la fonction 
définie par 



si a = , 

si a ^ est un carré de Z/pZ , 
si a n'est pas un carré dans Z/pZ 



Posons Y>(a) = et considérons la matrice entière symétrique M(n) 
avec coefficients 



Mij = V( 



i+j 



(mod p)) , < i, j < n . 
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Les propriétés du symbole de Legendre et (Eq. 1) montrent que M est 
p — autosimilaire . 

Exemple 2.1. Pour p = 2 la matrice M(oo) n'est rien d'autre que la 
réduction modulo 2 du triangle de Pascal (sous forme de matrice symétrique) 
avec coefficients dans {0, 1}. On a la factorisation 

*-Gi)-(ïï)(i-ï)00 

de la matrice de définition M et le Théorème 1.4 montre alors qu'on a 

n— 1 __ n— 1 

det(M(n))= 'fi (-l)^ Kl = I](- 1 ) Sfe 

où la suite (définie récursivement par so = 0, s 2 k = s^, et S2k+i = 
1 — Sfc) est la fameuse suite de Thue-Morse comptant (mod 2) le nombre 
de digits non-nuls d'un entier binaire (cf. [AS]). Plus précisément, on montre 
facilement qu'on a 

det(M(2n)) = (-1)" et det(M(2n + 1)) = (-l) n+s " . 

Le fait que la matrice LD = ^ j ^ ^ (et donc L(oo)D(oo)) soit 

d'ordre 2 est équivalente à la définition récursive suivante de la suite de 
Thue-Morse: sq = et s n G {0, 1} tel que 

E((fc) (mod 2)) (-1)-* = 

pour tout n > 1 (ici, (Q) (mod 2)) € {0,1} et l'égalité est sur Z). Ceci 
peut se généraliser comme suit: Le vecteur fi = (1, —1, —1, 1, . . . , (— l) Sn , . . .) 
associé à la suite de Thue-Morse vérifie 




î.e. fj, 



1 L ' j 



où L est la matrice 2— autosimilaire triangulaire inférieure définie par L 

(1 

1 1 



1 L'assertion (iii) du Théorème 1.4 montre donc que L 1 est 



2— autosimilaire, défini par L 1 = ^ j 1 / ' ^ rem P^ a( * an * ^ c i"d essus 
par la matrice 6— autosimilaire triangulaire inversible définie par Ljj = 
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G), < i, j < b, d'inverse R avec coefficients Rij = < i, j < b, 



on a 



L J = | 



i.e. /j' 



t -R\ 



pour le vecteur fi = (/io,/ii, • • •) = (1, —1, • • • ,/i n = (— . . .) avec n = 
J2i u ib l - En particulier, si 6 est impair, on a simplement fi n = (—1)™ ce qui 
se traduit par l'égalité 



E 

fc=0 (mod 2) 




k=l 




2 (E"*)-i 



pour n = E > 1 et fc = X) avec 6 > 3 un entier impair. On 
peut évidemment également considérer d'autres matrices b— autosimilaires 
triangulaires inférieures A à la place de la matrice L discutée ci-dessus et 
calculer la suite /x associée (qui n'est rien d'autre que la suite des coordonnées 
de la première colonne de la matrice b— autosimilaire A -1 ). 

Revenons maintenant au premier p = 2 et à notre matrice symétrique 
2— autosimilaire M définie par Mjj = Ç^) (mod 2) G {0, 1}. On peut 
montrer assez facilement que l'inverse M(n)~ 1 de la matrice M(n) (pour 
n G N) n'a que des coefficients dans {—1,0, 1}. 

Signalons finalement que les déterminants des sous-matrices de M(oo) 
obtenu en prenant l lignes et colonnes consécutives (et ayant donc comme 
coefficients Mj +s j +t , < i,j < l pour s,t € N fixé) ne semblent prendre 
que les valeurs ou ±1. 

Exemple 2.2. Pour le premier p = 3, la matrice M(oo) est la réduction 
du triangle de Pascal (mod 3) à valeurs dans {0, ±1}. La factorisation 
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montre que le déterminant 



n-l 



det(M(n))= [] W d ^i) 

(avec d(0) = 1, d(l) = (—2) et d(2) = — |) est une puissance de (—2). 
Example 2.3. Pour p = 5 on a M = LD (l) avec 
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et D une matrice diagonale avec coefficients diagonaux 

d(0) = 1, d(l) = -2, d(2) = 2, d(3) = -| d(4) = ^ . 
Le déterminant 

ra-1 

det(M(n)) = n n 

est donc un entier de la forme ±2 a 3 /3 avec a = a(ri) et (3 = (3 (ri) des entiers 
naturels convenables. 

Example 2.4. Pour le premier p = 7 la matrice M(oo) n'est plus non- 
dégénérée. On peut la perturber en considérant la matrice p— autosimilaire 
définie par 

1 

-1 
-1 
-1 
1 





M 



( 1 1 

1 t + 1 

1 -1 

f 1 

1 -1 

1 -1 

f 



1 
1 

-1 
-1 






1 1 1 \ 

-1 -f 

1 







o o o y 

qui est maintenant non-dégérée sur Z[t] et qui redonne la matrice initiale en 
posant t = 0. 

La matrice triangulaire L est alors donnée par 

/ 
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et la matrice diagonale D correspondante est donnée par 



d(0) = 1, d(l) = t, d(2) 
t - 10 



-2t-4 



, d(3) 



d(4) 



-, d(5) = 



-, d(6) = 



t + 2' 
-1 



4 ' ~ v ~' 2t-20 w t + 8 

Le déterminant det (M (ri)) se calcule maintenant en évaluant la fraction 
rationelle (qui est en fait toujours un élément de Z[t]) 



det (M (ri)) 



n—l 

n 



n 
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en t = 0. Il est soit soit de la forme ±2 a On peut montrer qu'il 
est nul si et seulement si m = H-fl 1 = n ~ 1 contient le chiffre \ii = 1 
ou si m contient le chiffre /ij = 2 suivi d'autre chose que de 666 . . . (i.e. 
det(M(l + J2 fJ-fl 1 )) = <^=^ ou bien il existe i avec [ii = 1 ou bien il existe 
i > j avec = 2 et fj,j ^ 6). 

Je remercie Jean-Paul Allouche pour ses commentaires et son intérêt. 

[AS] J-P. Allouche, J.Shallit, The ubiquitous Prouhet-Thue- Morse sé- 
quence, Proceedings of SETA 98 (C.Ding, T.Helleseth, H.Niederreiter, edi- 
tors), Springer (1999). 
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